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ABSTRACT 

As part of the TANAMl multi wavelength progam, we discuss new X-ray observations of the y-ray and radio-loud narrow line Seyfert 1 
galaxy (y-NLSl) PKS 2004-447. The active galaxy is a member of a small sample of radio-loud NLSls detected in y-rays by the 
Fermi Large Area Telescope. It stands out for being the radio-loudest and the only southern-hemisphere source in this sample. We 
present results from our X-ray monitoring program comprised of Swift snapshot observations from 2012 through 2014 and two new 
X-ray observations with XMM-Newton in 2012. Supplemented by archival data from 2004 and 2011, our data set allows for a careful 
analysis of the X-ray spectrum and variability of this peculiar source. The (0.5-10) keV spectrum is described well by a power law 
(T ~ 1.6), which can be interpreted as non-thermal emission from a relativistic jet. The source exhibits moderate flux variability on 
timescales of both months and years. Correlated brightness variations in the (0.5-2) keV and (2-10) keV bands are explained by a 
single variable spectral component, such as the one from the jet. A possible soft excess seen in the data from 2004 cannot be confirmed 
by the new XMM-Newton observations taken during low-flux states. Any contribution to the total flux in 2004 is less than 20% of 
the power-law component. The (0.5-10) keV luminosities of PKS 2004-447 are in the range of (0.5-2.7)x lO'*'* erg s“'. A comparison 
of the X-ray properties among the known y-NLS 1 galaxies shows that in four out of five cases the X-ray spectrum is dominated by 
a flat power law without intrinsic absorption. These objects are moderately variable in their brightness, while spectral variability is 
observed in at least two sources. The major difference across the X-ray spectra of y-NLS Is is the luminosity, which spans a range of 
almost two orders of magnitude from 10"*'* erg s“' to 10^ erg in the (0.5-10) keV band. 

Key words. galaxies:active - galaxies:individual(PKS 2004-447) - galaxies:jets - galaxies:Seyferts - quasars:general - X- 
rays:galaxies 


1. Introduction 

The recent detection of variable y-ray emission from radio-loud 
narrow line Seyfert 1 galaxies (y-NLS 1 hereafter) suggests that 
there are relativistic jets in these sources, which are similar to 


percent of the Eddington rate (see, e.g., |Marscher|2009 for a re¬ 
view). At the intersection of two rare types of AGN (NLSl and 
radio-loud AGN), y-NLS Is are unique laboratories for studying 
the jet formation in AGN. 


other types of active galactic nuclei (AGN), i.e. blazars and radio Eormally, NLS1 galaxies are classified by their strong and nar- 

galaxies (|Abdo et al.|2009b 

Eoschini|2012a[|D’Ammando et al.| row optical Balmer emission lines from the broad line region. 

2013bl Foschini et al. 2015 

1 . Typical properties of radio-quiet The intensities of the narrow lines are strong with respect to the 


(see, e.g.,|Hayashida 

2000( 

Wandel|2000(|Boroson|2002|l, an ac- 

cretion rat 
(see, e.g.. 

e close to t 
Foschini|2 

le Edc 
0T2b|l. 

ington limit, and spiral host galaxies 
Their detection in y-rays poses seri- 


ous questions as to the well-known paradigm that such power¬ 
ful jets are hosted by AGN in elliptical galaxies where a super- 
massive black hole (SMBH) of at least IO^Mq accretes at a few 


and [Oni]/Hy6 < 3 ( |Osterbrock & Pog^|1985| l. These prop¬ 
erties are typically accompanied by a strong Ee n emission line 
complex ( |Grupe|2004[ ) and strong characteristic high-ionization 
line emission from the broad line region. The origins of the 
strong Lett emission and of the narrow Balmer lines are still a 
topic of debate. In X-rays, compared to other Seyfert 1 galaxies. 
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NLS1 exhibit typically stronger X-ray flux variability on short 
timescales (see, e.g., Vaughan & Fabian 2003 [ Grupe 2004] Kara 


et al. 


1999 


Grupe|2004| l. 


2013 |l and a strong soft excess below 2 keV (Vaughan et al. 


Among NLS Is, only 7% are radio-loud (Komossa et al.pOOh 


Yuan et al. |2008|l, of which only a small number have 


been detected in y-rays ( Abdo et al.|POOP'S] |Foschini||2012a[ 
D’Ammando et ah] 2013b|l. Multi wavelength campaigns have 


been performed on three of the known y-NLSls detected by 
the Fermi Large Area Telescope (LAT; | Atwood et al.||2009l: 
PMNJ0948-1-0022 (|Abdo et al.||2009b|l, IH 0323-1-342 (|Paliya 
let al.|2014| l, and SBS 0846-1-513 ( [D’Ammando et al.|2013b| l7m 


present results from a multiwavelength monitoring of the radio- 
loudest source from this elusive y-NLSl sample. 

The active galaxy PKS 2004-447 is a y-NLSl galaxy at a 
redshift of z = 0.24 (|Drinkwater et al.| |1997|) with a ra- 


dio loudness parameter R = / 4 . 9 ghz//b ^ 1700 (Keller- 

mann et al. 

[19891 [Oshlack et al.[[2001[ hereafter OshOl; [Gallo 

et al. 

2006 

hereafter G06). Its optical properties are consis- 


tent with the NLSl classification: a flux ratio [Oin]/FIy8 =1.6 
and FWHM(H;0) = 1447 km s ’, even though the Feii emission 
{EWpea < lOA) is rather weak (OshOl). Based on the empirical 


relation between BLR size and optical luminosity (Kaspi et al. 
|2000| l, OshOl estimated the black hole mass from the H/3 line to 
be 10® ^Mq, using high-resolution data from the Siding Spring 
2.3 m Telescope. Studies in the X-rays (see e.g., G06; Paliya 
|et al.|20l'3] l suggest an unobscured spectrum in the (0.3-10) keV 
band, a possible soft excess below 1 keV, and mild long-term flux 
variability. 

We observed this source as part of the TANAMI multiwave¬ 
length monitoring program ( Ojha et al.|2010[ Kadler et al. 2015| l. 
Here we present results from the analysis of the X-ray spectra. 
In addition to archival data from 2004, we use two new XMM- 
Newton pointing observations performed in 2012 to explore the 
source’s long-term X-ray spectral variability. For the study of 
X-ray flux variability, PKS 2004-447 has been monitored with 
Swift since 2011. 

In a second paper (Schulz et al. submitted. Paper II), we re¬ 
port on the radio properties of PKS 2004-447, including the 
first TANAMI VLBI image at 8.4 GHz. VLBI imaging shows 
an elongated, one-sided radio jet extending from a bright com¬ 
pact core. Simultaneous multifrequency observations with the 
Australia Telescope Compact Array (ATCA) over four years re¬ 
veal a persistent steep spectrum. The estimated large-scale size 
of the radio emission and the spectrum are consistent with com¬ 
pact steep spectrum (CSS) sources, which are generally consid¬ 
ered to be young radio galaxies ( |0’Dea|1998| l. 

The paper is organized as follows. In Sect.|^ the monitoring pro¬ 
gram and data reduction methods are described. The analysis of 
X-ray data and its results are reported in Sect.|^and discussed in 
Sect.|^ which also includes a comparison with X-ray properties 
of other y-NLSls, CSS galaxies, and blazars. The main results 
of the work are summarized in Sect.|5] 


2. Observations and data reduction 

PKS 2004-447 has been monitored in X-rays since May 2013 
as part of the TANAMI multiwavelength program. We used the 
space observatories XMM-Newton pansen et al.|2001]l and Swift 


(Gehrels 


2004 |l, which provide X-ray data and optical/UV cover¬ 


age. High-angular-resolution radio observations have been per¬ 
formed since October 2010 by the TANAMI VLBI program. In 
addition, the source has been observed several times by ATCA. 
Simultaneous y-ray data are provided by the FermifLAI. The 
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X-ray monitoring consists of two XMM-Newton pointings in 
2012 separated by five months and connected by four Swift ob¬ 
servations. We have continued observations with Swift in 2013 
and 2014. To study the long-term behavior of PKS 2004-447, 
archival data from XMM-Newton and Swift prior to 2013 May 
have been included in the analysis. Detailed information on all 
observations is provided in Table 


2.1. XMM-Newton observations 


We observed PKS 2004-447 twice with XMM-Newton in 2012: 
for 37.9 ks on 2012 May 1 and for 39.9 ks on 2012 October 
18, hereafter referred to as X-2012-05-01 and X-2012-10-18. 
The archival XMM-Newton observation from 2004 April 11 (for 
41.9 ks, X-2004-04-11 hereafter) has been discussed in detail by 
G06. All observations were performed with the European Pho¬ 
ton Imaging Camera (EPIC) using PN ( IStriider et al.|200T] ) and 
MOS ( |Turner et al.|2001a I CCD arrays and the Reflecting Grat¬ 
ing Spectrometer (RGS, den Herder et al.|2001| l, as well as the 
Optical Monitor (OM, jMason et al.||2001| . The EPIC operated 
in full-frame mode with the medium filter in 2004 and the thin 
filter in 2012. Observation data files (ODEs) were processed to 
create calibrated event lists and full frame images using standard 
methods with the XMM-Newton Science Analysis System (SAS 
Version 12). Lightcurves were extracted for the (0.5 - 2)keV, 
(2 - 10)keV, and (0.5 - 10)keV energy ranges, adopting a cir¬ 
cular extraction region with 30 arcsec radius for the source. The 
background regions were chosen from a source-free circular re¬ 
gion of the same radius on the same chip with similar distances 
to the readout nodes to ensure similar background noise levels. 
Screening for flaring particle background was done following 
the method to maximize the signal-to-noise ratio (S/N) by |Picon-] 
celli et al. ( 2004| l. Pile-up was negligible in all observations. The 
net exposure times are listed in Table [T] Source and background 
spectra were extracted for single and double event patterns from 
the filtered event lists. Response files were generated using rmf- 
GEN and ARFGEN of the SAS software package. The RGS data had 
insufficient counts for a dedicated spectral analysis. 


2.2. Swift observations 

We monitored PKS 2004-447 with the Swift X-ray telescope 
(XRT, [Burrows et al.|2005| ) from the beginning of 2012 through 
the beginning of 2014. In addition, we used archival data from 
2011, which were discussed by [Paliya et al.| ( [2013| l. Details 
of each observation are listed in Table [IJ Three observations 
had relatively short exposures (2.3 ks, 1.6 ks, and 0.6 ks). Ow¬ 
ing to the low S/N in the respective X-ray spectra, they are not 
considered in our analysis. Eor all other observations, we used 
data from the XRT, which operated in photon-counting mode 
(PC). The data were cleaned, and calibrated event files were 
created using the standard filtering methods and the xrtpipeline 
task, distributed in the HEASARC (v6.12.0) within the HEASoft 
package. Source and background spectra and lightcurves in the 
(0.5 - 10) keV energy range were created, using a circular source 
region with a radius of 35 arcsec. The net exposure times are also 
listed in Table [T] To increase the S/N of the spectra. Swift obser¬ 
vations that were observed within a few days of each other were 
merged after checking that the spectra did not show significant 
flux or spectral variability between the averaged epochs. 








































































































A. Kreikenbohm et al.: The gamma-ray emitting radio-loud narrow-line Seyfert 1 galaxy PKS 2004-447 


Table 1. Details of all XMM-Newton and Swift observations of PKS 2004-447. 


Inst-ObsDate 

Obsid 

Duration 

Net exposure 

Spectrum counts 



)ks] 

)ks] 

)counts] 

(1) 

(2) 

(3) 

(4) 

(5) 

X-2012-10-18 

0694530201 

39.8 

34.5 

10788 

X-2012-05-01 

0694530101 

37.9 

36.5 

6125 

X-2004-04-11 

0200360201 

41.9 

38.4 

15831 

S-2014-03-14 

000324920) 16,17]“ 

17.2 

17.1 

248 

S-2013-11-19 

000324920)14,15]* 

17.7 

16.7 

523 

S-2013-10-13 

0003249200)09-13]“ 

16.5 

16.5 

302 

S-2013-09-27 

00032492007 

8.3 

8.3 

234 

S-2013-07-07 

0003249200)5,6]'^ 

23.1 

23.0 

444 

S-2012-09-30 

00032492004 

6.0 

6.0 

63 

S-2012-09-12 

00032492003 

5.4 

5.4 

63 

S-2012-07-22 

00032492002“ 

2.3 

- 

- 

S-2012-07-03 

00032492001 

4.9 

4.9 

67 

S-2012-03-14 

00091031007 

7.3 

7.3 

60 

S-2011-11-15 

00091031006 

7.1 

7.1 

73 

S-2011-09-17 

00091031005 

7.5 

7.34 

167 

S-2011-07-29 

00091031004“ 

1.6 

- 

- 

S-2011-07-25 

00091031003“ 

0.6 

- 

- 

S-2011-07-14 

00091031002 

5.0 

4.9 

56 

S-2011-05-15 

00091031001 

6.8 

6.8 

61 


Notes. (1) Instrument and observation date. We use the format Inst-yyyy-mm-dd: X - XMM-Newton EPIC, S - Swift XRT. (2) Observation 
identifier. (3) Observation duration. (4) Net exposure time after screening for flaring particle background. (5) Background-subtracted total num¬ 
ber of counts in the (0.3-10) keV band. Merged observations: 00032492016 (2014-03-14, 7.6 ks) and 00032492017 (2014-03-16, 9.7 ks). 
® 00032492014 (2013-11-19, 5.5ks) and 00032492015 (2013-11-20, 12.2ks). 00032492009 (2013-10-13, 4.7ks), 00032492010 (2013-10- 

20, 4.11 ks), 00032492011 (2013-10-27, 2.35 ks), 00032492012 (2013-10-29, 1.51 ks), and 00032492013 (2013-11-03, 3.88ks). ^ 00032492005 
(2013-07-07, 11.4ks) and 00032492006 (2013-07-14, 11.6ks). Because of the relatively short exposure time, this observation was not included 
in the analysis. 


3. X-ray analysis 

Spectral fitting was done using the Interactive Spectral In¬ 
terpretation System (ISIS, Version 1.6.2-30, |Houck & Deni-| 
cola 2000) 1. Unless stated otherwise, uncertainties correspond 
to 90% confidence limits on one parameter of interest 
- 2.7). In the following, the standard cosmological model 
{Hq - 70kms“* Mpc“\ Qm = 0.3, and A = 0.7) is assumed. 


3.1. XMM-Newton PN/MOS data 


Visual inspection of the XMM-Newton data revealed a flat and 
smooth spectrum for each observation (see Figs. [T]and[^. The 
source is detected from 0.5 keV to lOkeV. We first adopted 
an unbroken power law modulated by Galactic HI absorption 
Ah.GbI = 3.17 X 10^° cm“^, based on the Galactic Hydrogen LAB 
survey ( jKalberla e t al.|2005]l. The high- resolution ISM absorp¬ 
tion model tbnev\{^ ( Wilms et al. 2000 1 is used to account for 
neutral absorption. Cross-sections and abundances were taken 
from |Verner et al.| ( | 1996) 1 and )Wilms et al.) ( 2000| l, respectively. 
The large number of spectral data counts allowed for the use of 
;^^^-statistics, and each spectrum was grouped to a S/N of five per 
energy bin. For each observation, we fit the available PN and 
MOS data in the (0.5-10) keV energy range simultaneously. 

The absorbed power-law model yielded a good fit to each spec¬ 
trum. Best-fit results for the photon index F lay between 1.50 and 
1.65. We found no evidence of intrinsic absorption when leav¬ 
ing the column density of the absorption free to vary. The upper 
limit of the column density was Ah < L7AH,Gai- F was there- 


‘ Wilms, J., Jutt, A. M., Schulz, N. S., Nowak, M. A. (2012), published 
at http://pulsar.sternwarte.uni-erlangen.de/wilms/research/tbabs/ 


fore fixed to Ah.GqI to reduce the number of free parameters. In 
the following, unless stated otherwise, all spectral models men¬ 
tioned in this paper include the Galactic absorption component 
with fixed column density. The XMM-Newton spectra are shown 
in Figs.[^and[^ along with their best-fit models (solid line). We 
convolved the power-law component of the best-fit model with 
cflux to determine the unabsorbed energy fiux in different energy 
bands, i.e. total band (0.5-10) keV, soft band (0.5-2) keV, and 
hard band (2-10)keV. For each spectrum, we performed three 
independent fits by taking only data within each respective band 
into account. The sub-bands were defined based on typical fiux 
ranges used in the literature to facilitate comparison. The best-fit 
parameters and absorption-corrected fiuxes of the fit are listed in 
Table ]2] 

We found no indications of any spectral line emission features. 
The addition of a Gaussian line at 6.4keV with a fixed width 
of 1 eV did not yield any significant improvement of the fit 
(Ax^ = -0.1 for one additional free parameter). We measured 
an upper limit for the equivalent width of FWe^kev ^ 60 eV in 
all three observations. This is consistent with the results of G06. 
G06 also report indications of a weak soft excess below ~ 1 keV. 
They found that a broken power law with a break energy at 
~ 0.6 keV is a good approximation of the data. Thus, we looked 
for the existence of such a soft excess by modeling the spectra 
only above 2keV and extrapolating the best fit down to 0.5 keV. 
No significant excess emission is observed in the case of X- 
2012-05-01 or X-2012-10-18, but visual inspection of the resid¬ 
uals suggests a possible soft excess in the data of X-2004-04-11 
(see Figs. [^and]^. The spectra of all three XMM-Newton obser¬ 
vations were then modeled phenomenologically with a broken 
power law. We note that results for the soft photon index Fi are 
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Table 2. Best-fit results of the absorbed power-law model for XMM-Newton and Swift spectra. 


Inst-ObsDate 

( 1 ) 

X-2012-10-18 

X-2012-05-01 

X-2004-04-11 

S-2014-03-14 

S-2013-11-19 

S-2013-10-13 

S-2013-09-27 

S-2013-07-07 

S-2012-09-30 

S-2012-09-12 

S-2012-07-03 

S-2012-03-14 

S-2011-11-15 

S-2011-09-17 

S-2011-07-14 

S-2011-05-15 

S-low 

S-medium 

S-high 


ObsID 

( 2 ) 

0694530201 

0694530101 

0200360201 

000324920[16,17]' 

000324920[14,15]'' 

000324920[09-13]“ 

00032492007* 

0003249200[5,6]'" 

00032492004' 

00032492003' 

00032492001' 

00091031007' 

00091031006' 

00091031005“ 

00091031002' 

00091031001' 


r 

(3) 


^■^^-0.02 


).05 


1.56_ 

1 52“^^-^^ 

^■^■^-0.12 
1 64+"'^ 
52^8:1? 
-Q.13 
1 47+0-17 

1.41 
1.63 
1.77; 
1.75 
1.57; 
2.04; 
1.45 
1.66 
1.38 


■fO.12 




3 


1 53+0 

1 50+"'° 

1 57+°'^ 

+ -n.i 2 


-P’o.S-lOkeV 

(4) 


0 70 + 0-02 

^-Q-Q2 

0.48+° 

1 45+° 

1.52+; 
0.99^; 
1.54+; 
1.16+; 
0.52+1 
0.56i; 
0.63+; 
0.42+1 
0.39^; 

0.63+; 
0.54+: 


0.53!3 

1-02!| 




-P’o.5-2keV 

(5) 


7y23+injr 


0.15+' 

0.43 


01 

8 : 8 ? 

01 


0.26^; 

0.51!; 

0.30!; 

0.43!; 

0.31!; 

0.16!; 

0.22!; 

0.24!; 

0.13!; 

0.18!; 

0.33!; 

0.21!; 

0.14+; 


oi^ 

0.30!; 

o-46!8 


-P’2-10keV 

(6) 


04^ 

0.35!|»^ 

1.06^ 


;8:i^ 

;8:f8 


0.75! 
0.91! 
0.71 
1.21 
0.92 
0 26^«:‘S 

-Q-15 

0 41+0.18 
'+■^ + -0.30 

0 50+°-^' 
- 0-35 

0 31+° '^ 
-Q19 

0 19+0.08 
'+■ + ^-0.14 
0 79+0-21 
'^-0.28 
0 49+0-18 

0 56+°-^' 

^•~”+-n.34 

0 33+° °' 

- 0-12 
0 70+0-09 
-0.14 
0 94+0-15 


7^0.5-lOkeV 

(7) 

1.11+0.03 
0.75 ± 0.04 
2.26 + 0.05 
1.45 + 0.23 
2.44 + 0.26 
1.55 + 0.23 
2.38 + 0.42 
1.77 + 0.22 
0.82 + 0.25 
0.93 + 0.26 
1.03 + 0.29 
0.67 + 0.21 
0.68 + 0.17 
1.82 + 0.37 
1.02 + 0.32 
0.82 + 0.27 
0.84 + 0.11 
1.59 + 0.16 
2.36 + 0.28 


stat/dof (statistic) 

( 8 ) 

296.4/258 ix^) 
135.3/131 ix^) 
388.3/398 (x^) 
622.1/948 (C) 
594.7/948 (C) 
517.2/948 (C) 
481.8/948 (C) 
671.6/948 (C) 
244.8/948 (C) 
270.2/948 (C) 
249.8/948 (C) 
271.9/948 (C) 
243.5/948 (C) 
405.2/948 (C) 
258.2/948 (C) 
273.2/948 (C) 
15.7/12 ix^) 
19.7/27 0^2) 
22.6/23 (x^) 


Notes. Best-fit parameters for the absorbed power law. The absorption is fixed to its Galactic absorption, 3.17 x 10^°cm“^, based on the survey 
by |Kalberla et al. ( |2005[ l. (1) Instrument and observation dates as in Table fTj (2) Observation ID, [n,m] indicates merged spectra. (3) Power-law 
photon index. (4,5,6) Absorption-corrected fluxes in units of 10“*^ erg cm“^“‘ for the (0.5 - 10) keV, (0.5 - 2) keV, (2 - 10) keV energy ranges. 


Swift data are not rebinned. Since the degrees of freedom are then given hy the energy grid of the XRT it is constant for all Swift observations. 
(i.m.h) Observations included in the merged Swift data set of low (1), medium (m), or high (h) fluxes (see Sects. lTn & l3.2l for details). 


not well constrained because of the narrow energy range from 
0.5 keV to the break energy. Fit results for the broken power law 
are summarized in Table For X-2004-04-11, the model pro¬ 
vided a slightly better fit than the simple power law. The best- 
fit parameters are consistent with results by G06, except for the 
soft photon index Fi, which cannot be constrained well. As for 
X-2012-05-01, it was not possible to constrain either of the pho¬ 
ton indices when we left the break energy free to vary. We fixed 
the break energy to 0.6 keV as reported in G06, which allowed 
us to obtain constraints on Fi and F 2 . The fit of X-2012-10- 
18 yields good constraints, however the break energy is signifi¬ 
cantly higher than in G06. Within its uncertainties, the hard pho¬ 
ton index F 2 is roughly consistent with the one from the simple 
power law in all three observations. 

We tested the significance of the broken power-law model by 
applying the F test. The null-hypothesis model assumes an un¬ 
broken power law as an underlying model. A broken power-law 
model was considered to be significantly different from the sim¬ 
ple power-law model if the probability that the latter is correct 
was less than 1% (expressed by the p value). The F-test p values 
are also listed in Tabled For X-2012-05-01 and X-2012-10-18, 
both models were statistically equal. For X-2004-04-11, the bro¬ 
ken power law gave a slightly better fit than the simple power 
law, and the p-value was inconsistent with the null hypothesis 
according to our threshold. Thus, the presence of a weak soft 
excess is tentative. However, it was striking that no excess was 
detected in the observations in 2012. Since background count 
rates were consistent with each other in all three observations, a 
contribution of systematic background variations to the soft X- 
rays could be excluded. 
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Assuming the soft excess of X-2004-04-11 is caused by an¬ 
other emission component, the spectrum was modeled phe¬ 
nomenologically using a blackbody in addition to the simple 
power law. The power law was fitted in the (2-10) keV band 
and subsequently fixed to the best-fit parameters (see Table 
1^. This model yielded a good fit of ;^^^/d.o.f. = 396.4/398. 
The contribution of the blackbody to the (0.5-2) keV flux is 
(0.06 + 0.01) X 10^'^ ergcm^^ s“', which makes up 10-20% of 
the power-law flux in the soft band. We applied the same model 
to the data of 2012 (see Table where no soft excess has been 
observed. The fits were acceptable and yielded a contribution 
of the blackbody to the soft band between 5% and20% of the 
power-law flux in the same band. We note that for these fits, the 
power law is slightly steeper, which may in principle hide the 
presence of a faint extra component. 

Two difference spectra were created by subtracting the 
background-subtracted EPIC PN spectra of X-2012-10-18 and 
X-2012-05-01 from the one of X-2004-04-11 (see Fig.j^. The 
spectra were rebinned in the same way as the original data. We 
ignored energy bins larger than 1 keV. This leads to no spec¬ 
tral bins being regarded above 7keV, because the number of 
spectral data counts of the X-2012-05-01 spectrum decreases 
rapidly in this range. We fit the difference spectra by apply¬ 
ing a Galactically-absorbed simple power law to each of the 
(0.5-7) keV and (2-7) keV bands independently (see Table [^. 
In the total range, the best-fit parameters of both difference spec¬ 
tra are consistent with those of X-2012-10-18 and X-2012-05- 
01 within 90% uncertainties. The best-fit values of the photon 
indices are slightly flatter. This is expected if the spectrum of 
X-2004-04-11 is harder than in 2012. Fitting only the hard band 
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Energy [keV] 



Fig. 1. XMM-Newton spectrum and best-fit results of X-2004-04-11 
(top) and X-2012-05-01 (bottom). Upper panel a): EPIC pn (blue) 
and MOS (gray) data, as well as corresponding best fits for an ab¬ 
sorbed power law evaluated over (0.5-10) keV (black solid line) and 
(2-10) keV (black dashed line). Mid panel b): Data-to-model ratio for 
the best fit over the whole (0.5-10) keV energy range. Lower panel c): 
Respective residuals for fitting the (2-10) keV energy range alone, ex¬ 
trapolated to lower energies. A tentative soft excess below 0.6 keV is 
only observed for X-2004-04-11 (compare shaded regions). 


and extrapolating the fit into the soft band, a soft excess is not 
observed in the difference spectra, as illustrated in Fig. 


3.2. Swift/XRT 

The Swift XRT spectra have significantly lower S/N, and in 
most observations the source is not detected above 6keV. Be¬ 
cause of the poor quality of the data, we did not adopt a 



; _ ^ ^^ ^ ^ ^^^^^^^^ 

0.5 r 2 5 ro 

Energy [keV] 


Fig. 2. XMM-Newton spectrum and best-fit results of X-2012-10-18. 
Upper panel a): EPIC pn (blue) and MOS (gray) data, as well as cor¬ 
responding best fits for an absorbed power law evaluated over (0.5- 
10)keV (black solid line) and (2-10) keV (black dashed line). Mid 
panel b): Data-to-model ratio for the best fit over the whole (0.5- 
10)keV energy range. Lower panel c): Respective residuals for fitting 
the (2-10) keV energy range alone, extrapolated to lower energies. A 
tentative soft excess below 0.6 keV is only observed for X-2004-04-11 
(compare shaded regions). 



L_^^^^ _I_^^^^^^^^ 

0.5 1 2 5 10 

Energy [keV] 


Fig. 3. XMM-Newton ratio of the XMM-Newton EPIC PN difference 
spectra to a power-law model. Difference spectra were created by sub¬ 
tracting the background-subtracted spectra of X-2012-05-01 (blue) and 
X-2012-10-18 (black) from the one of X-2004-04-11. Simple power- 
law fits were performed on the energy ranges (0.5-7) keV and (2.0- 
7) keV, where the fit to the latter band was then extrapolated to lower 
energies. 
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Table 5. Power-law fit of the XMM-Newton EPIC PN difference spectra. 


Subtracted spectrum 

r2-10keV 

TO-S-lOkeV 

T"/d.o.f. 

X-2012-10-18 

1.36 + 0.26 

1.43+0.09 

29.4/33 

X-2012-05-01 

1.23+0.18 

1.45+0.06 

103.6/69 


Notes. Best-fit parameters for the simple power-law fit in the 
(0.5-10) keV and (2.0-10) keV energy hand to the difference spectra. A 
difference spectrum is computed hy subtracting the respective spectrum 
from X-2004-04-11. 


more complex model than the simple power law, which al¬ 
ready gave a good description of the XMM-Newton data. The 
model is fitted to the available unbinned data in the (0.5- 
10)keV range using the unbinned likelihood statistic Cstat 
( |Cash| [T979| ). Best-fit parameters are listed in Table The 
low number of spectral data counts mean that the fit parame¬ 
ters were affected by high uncertainties, especially during times 
of low flux. To improve statistics of the spectral fit, we de¬ 
fined three flux states based on the results of the single ob¬ 
servations and merged spectra with fluxes within these inter¬ 
vals (in units of 10“'^ erg cm“ 2 s ^); 0.4 < Fo.s-ikev ^ 0.9 (low), 
0.9 < T'o. 5 - 2 keV < 1-4 (medium), and 1.4 < T’o. 5 - 2 keV (high). This 
gave a total exposure of 57.02 ks with 837 counts in the low- 
flux state, 46.9 ks with 996 counts in the medium-flux state, and 
25.1 Iks with 795 counts in the high-flux state. The quality of 
the merged spectra then enabled the use of ;^f^-statistics and the 
spectra were analyzed analogously to the XMM-Newton spectra. 
For each of the three spectra, simple power-law models yield 
good fits with no significant evidence of any other feature such 
as a soft excess. Best-fit parameters for these fits are listed in 
Table]^ labeled as S-low, S-medium, and S-high. There was no 
evidence of flux-dependent spectral variations. 


3.3. Multitimescale X-ray variability 


In the following, reported fluxes and luminosities refer to their 
absorption-corrected values. If not stated otherwise, source 
count rates refer to their background-subtracted values. Figure 
|4]shows the evolution of the total flux of the source over time. 
PKS 2004-447 shows moderate X-ray flux variability on 
timescales of weeks to years. The scheduling of the monitor¬ 
ing observations does not allow for detecting variability on 
timescales shorter than two weeks. For the XMM-Newton ob¬ 
servations from 2004, the total flux is about 50% lower in those 
from 2012 and regained a similar value in September 2013. 

We calculated the fractional variability amplitude Avar following 
the definition by [Vaughan et al.| ( |2003| l: 


A 


var 



err (Avar) 



( 1 ) 

( 2 ) 


where N is the number of data points in the time series with a 
mean of the fluxes x. The variance and mean square error of the 

time series are denoted as 5 ^ and crg^r, respectively. 

For the full (0.5-10) keV energy range the fractional variabil¬ 
ity is A™'*' = 0.44 + 0.19. If divided into the soft (0.5-2) keV 
and hard (2-10) keV bands, we find A*^'= 0.39 + 0.18 and 
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Av“‘* = 0.35 + 0.33, consistent with the fractional variability of 
the total range. As a result, both the soft and hard ranges show 
the same variability characteristics. 

For a model-independent analysis of spectral variability, we de¬ 
termined the number of background-subtracted source counts in 
the soft and hard bands. Since source counts are not affected by 
calibrations issues at low energies, we extended the soft band 
to (0.2-2.0)keV. One-sigma uncertainties, ctn, of the number of 
counts, N, in a given band were approximated followingjGehrels 
( [T986l l: 

o-Ar=s VV-hO.75 h- 1. (3) 

The number of detected counts strongly depends on the 
energy-dependent effective area of a given instrument. Thus, 
for similar fluxes we expect different count rates and a number 
of spectral counts for Swift XRT and XMM-Newton EPIC PN. 
There were no simultaneous observations of both Swift and 
XMM-Newton, which could have been used to calculate the 
conversion factor between these instruments. Thus, counts were 
analyzed separately for each telescope. 


Hardness ratios HR were derived using the fractional difference. 


S -H 

HR = 

S +H 

(4) 

cthr = 2(S +H)-^ ^H^crj+S^crl 

(5) 


where S and H are source counts in the soft and hard bands with 
one-sigma uncertainties, cts and cth, respectively. 

Figure]^ shows the hardness ratio plotted as a function of the to¬ 
tal source count rate (hardness intensity diagram, HID) for the 
different instruments. For both instruments, the data are consis¬ 
tent with a constant HR. 

Spectral analysis of the XMM-Newton data yielded hints of a 
second spectral component visible below 2keV. We tested the 
presence of a non- or a less-variable component by comparing 
the variation in the count rates in the soft and hard bands (see 
Fig. 1^ for the hard versus soft band diagram, HSD). There is a 
wide spread at lower count rates. However, the data are overall 
consistent with a linear relation between the soft and hard bands. 
A linear dependency between these bands suggests that the flux 
variation of the spectrum is dominated by those of a single com¬ 
ponent. The presence of a less varying component would induce 
a positive or negative offset from the origin. We fitted both the 
whole data set and only high-S/N data with a linear regression 
and found a relation fully consistent with a line through the ori¬ 
gin (see Table 1^. The fits were consistent with a single compo¬ 
nent at 2 (T for all instruments. 

A comparison of background-subtracted count rates of the indi¬ 
vidual and merged spectra shows that spectra with lower counts 
appear softer. 


4. Discussion 

4.1. The X-ray spectrum of PKS2004-447 

The X-ray emission of PKS 2004-447 in the (0.5-10keV) en¬ 
ergy range is described best by a simple, Galactically-absorbed 
power law. A tentative soft excess is observed in the XMM- 
Newton spectrum of 2004. However, this excess is not detected in 
our new observations from 2012. The difference spectrum is de¬ 
scribed well by a simple power law, but the S/N is not sufficient 
to draw any conclusion about the tentative excess (see Fig.|^and 


12 






























A. Kreikenbohm et al.: The gamma-ray emitting radio-loud narrow-line Seyfert 1 galaxy PKS 2004-447 
Table 3. Results of the absorbed broken power-law model for XMM-Newton spectra. 


Inst-Obsdate 

Ti 

Tz 

Eb 

T^d.o.f. 

p-value 

(1) 

( 2 ) 

(3) 

(4) 

(5) 

(6) 

X-2012-10-18 

1 ^2 0 + 0.16 

1 

^•^'^- 0.01 

1 02 +°-^^ 
^•'^^- 0.16 

288.6/256 

0.033 

X-2012-05-01 

2 4Q + 1.60 

^^^- 0.06 

0.61 (fixed) 

134.5/130 

0.381 

X-2004-04-11 

> 3.05 

^^^- 0.03 

0 

^•'^^- 0.05 

372.1/396 

0.0002 


Notes. Best-fit parameters for the absorbed broken power law. The absorption is fixed to its Galactic value, 3.17 x 10^° cm“^ ( jKalberla et al. 2005 i. 
(1) Instrument and observation date. (2) Soft photon index below £b- (3) Hard photon index above Eb- (4) Break energy in units of keV. (5) Fit 
statistics of the broken power-law model. (6) F-test p value, i.e. likelihood that the null-hypothesis model is correct (see Sect. |3.1| for details). 


Table 4. Results of the black body and power-law model for XMM-Newton spectra. 


Inst-Obsdate 

Tz-lOkeV 

kT 

Ebb 

FpL 

X^/d.o.f. 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

X-2012-10-18 

1.57 + 0.10 

^■■^^- 0.12 

0.012 + 0.005 

0.216 + 0.005 

291.9/258 

X-2012-05-01 

1.43 + 0.19 

A 90 + 0.06 

^*^'^-0.07 

0.023 + 0.005 

0.128 + 0.004 

134.3/131 

X-2004-04-11 

1.39 + 0.08 

^*^'^-0.07 

0.059 + 0.009 

0.383 + 0.006 

396.4/398 


Notes. Best-fit parameters for the absorbed simple power law and black body. The absorption is fixed to its Galactic value, 3.17 x 10^°cm“^ 
( [Kalberla et al.|2005| l. (1) Instrument and observation date. (2) Photon index, fitted in the (2-10) keV band and fixed for broad band analysis. (3) 
Black body temperature in keV. (4,5) Absorption-corrected (0.5-2) keV flux in units of 10“*^ erg cm“^ s“' of the black body and the power law, 
respectively. (6) Fit statistics'^ per degrees of freedom. 


Year 

2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 



Fig. 4. Long-term behavior of the total absorption-corrected 0.5-10 keV flux. The gray-shaded background denotes the multiwavelength monitor¬ 
ing period. XMM-Newton data are shown as black squares. Swift data are plotted as triangles. 


Table 6. Results from the linear regression to the HSD. 


Data set 

Slope a 

Offset [counts s '] 

XMM-Newton 

2.26 + 0.04 

0.027 + 0.006 

Swift 

1.84 + 0.33 

0.004 + 0.002 

Swift merged 

2.10 + 0.43 

-0.001 + 0.002 


Notes. Best-fit parameters for the linear regression Csoft = rrChard + b, 
where Csoft and Chard are the spectral counts in the soft and hard bands, 
respectively. Uncertainties correspond to 1-sigma uncertainties. 


Table |5]l. The power-law spectrum and photon index are remi¬ 
niscent of the non-thermal X-ray emission of type-1 AGN. The 
spectrum exhibits flux variability across the entire band. Com¬ 
paring best-fit models of each observation, we find evidence of 
weak spectral variability, while there is no evidence of any in¬ 
trinsic absorption. This suggests that the power-law component 
itself varies, caused by changes in the emission region. Although 
we do see evidence of weak spectral variability from spectral 
fitting, the variations between the source brightness and spec¬ 
tral shape are not correlated (Fig. |^. The linear behavior of the 
count rates in the HSD (Fig.[^ can be explained by a dominating 
power-law component that exhibits changes in its normalization. 
These variations cause a homogeneous flux variation across the 
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Fig. 5. Hardness-intensity diagram (HID): Hardness ratio S-H/S+H as a function of total background subtracted source count rate. Error bars 
correspond to 1-sigma uncertainties. Left: Results from Swift observations. Individual observations are shown as gray open symbols, while filled 
symbols denote results from merged data sets. Right: HID of the XMM-Newton data. 
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Fig. 6. Left: Hard- versus soft-band diagram (HSD): Comparison of the count rates in the soft (0.2-2) keV and hard (2-10) keV band. Error bars 
correspond to 1-sigma uncertainties. Lines denote the best fit of the linear regression for low- (dashed) and high-S/N (solid) data. Left: Results 
from Swift observations. Individual observations are shown as gray open symbols, while filled symbols denote results from merged data sets. 
Right: HSD of the XMM-Newton data. Dashed line shows the best fit of the linear regression. 


(0.5-10) keV band. For interpreting the X-ray spectrum, we con¬ 
sider the following three scenarios: 

(I) a non-jet-dominated X-ray spectrum (characteristic of Seyfert 
galaxies), for which the unbroken power-law spectrum is asso¬ 
ciated with thermal inverse-Compton scattering of optical and 
UV photons from the accretion disk in the disk corona. A soft 
excess in these sources may be caused by relativistically blurred 
reflection from the accretion disk ([Crummy et al.|200^ or ther¬ 


mal Comptonization of disk emission by a population of elec¬ 
trons with low temperature and high optical depth (see, e.g., 
[Haardt & Marasclil |1993| ). We note that the photon index of 
PKS 2004-447 is significantly harder than those observed in 
radio-quiet NLS1 s (see, e.g.,[Vaughan et al.|1999[|Zhou & Zhang 
[20T0l[Grupe et al.|2010[ l. 


(II) a Seyfert-Iike X-ray spectrum that may be contaminated by 
a significant contribution from a relativistic jet. In this scenario, 
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jet-synchrotron emission contributes to the soft X-ray spectrum 
in 2004, while the unbroken power law can be associated with 
the Seyfert-like component. The detection in y-rays, however, 
suggests that the jet dominates at higher energies. Thus, we do 
not favor any scenario where the jet emission only contributes in 
the form of a soft excess. 

(Ill) a jet-dominated X-ray spectrum, characteristic ofblazars, in 
which the spectrum is associated with X-ray emission from the 
jet alone. Emission from the accretion disk is negligible. Within 
the framework of leptonic models, blazar X-ray spectra may be 
comprised of both synchrotron emission and inverse-Compton 
emission from relativistic particles in the jet plasma (see, e.g.. 


Maraschi et al. 

119921 [Dermer & Schlickeiser 

[I993l [Bloom & 

Marscher] 1996 

Ghi^llini et al.| 1998 2010 

[Bottcher et al. 


2013^. In this scenario, the hard power law can be interpreted as 


inverse-Compton emission from the jet and synchrotron emis¬ 
sion would account for the tentative soft excess in the data of 
2004. 

G06 have modeled the SED of PKS 2004-447 from the ra¬ 
dio band to hard X-ray band using a simple one-zone syn¬ 
chrotron self-Compton (SSC) model. They find that the X-ray 
spectrum could not be explained by a single SSC component, 
but an additional weaker emission component of (1.2 + 0.4) x 
10“^"^ ergs^'cm^^keV^' at 2keV is required. They argue that the 
X-ray spectrum is dominated by thermal Comptonization in ad¬ 
dition to a weak SSC component that gives rise to the tentative 
soft excess. |Paliya et al.| ( |2013[ ) also modeled the SED includ¬ 
ing y-ray data using a jet-emission model that consisted of an 
SSC and an external Compton component. A thermal Comp¬ 
tonization component was not needed. They further conclude 
that the SED of PKS 2004-447 resembles the distribution of flat- 
spectrum radio quasars (ESRQ). Eor these sources, hard X-ray 
emission above 2 keV is interpreted as inverse-Compton radia¬ 
tion, which suggests a non-thermal origin for the X-ray spectrum 
of PKS 2004-447. 


The jet-like non-thermal origin for the dominating X-ray emis¬ 
sion is supported by the flat photon index with respect to radio¬ 
quiet NLSls. This is consistent with values typically found in 
blazars, particularly ESRQs (see, e.g.. Pan et al. |2012[ [Rivers | 


jet al.|2013| [Dai & Zhang|2003[ jSambruna et al.|2004||. In Paper 
II we analyze the TANAMI VLBI data of PKS 2004-447 and 
And a one-sided core-jet radio structure. Together with the X-ray 
data, this result favors the third scenario and strongly implies that 
emission from the relativistic jet is dominating the X-ray spec¬ 
trum and that additional thermal or non-thermal emission from 
the accretion disk and its corona may only play a role in very 
low flux states of the source, as indicated by the HSD (see Pig. 
1 ^. 


4.2. Soft excess 

A faint soft excess is detected only tentatively in the XMM- 
Newton observation X-2004-04-11, consistent with findings 
by G06. This component could not be verified by our new 
XMM-Newton observations in 2012. Pitting a power law to 
the XMM-Newton data only in the (2-10)keV range (see 
dashed lines in Pigs. and leads to an excess emission of 
Esx = (6 + 1) X lO^^'^ergs^'cm^^ in the (0.5-2) keV band (see 
Sect.[3lTl. 

It remains questionable whether it is appropriate to interpret this 
as an additional spectral component. A variable soft component 
may explain the lack of excess emission in the data of 2012. 
We note that in radio-quiet Seyferts that include a few NLSls, 
there have been cases where the soft excess varies independently 


from and more slowly than the hard power law ( 

Turner et al. 

2001b Edelson et al.[[2002[ Mehdipour et al.[[201 

I| Arevalo & 

Markowitz 2014|l, and even seems to disappear 

fom one X- 


& Reeves 2009 Rivers et al.||2012l. Such a component should 


cause deviations from a linear function through the origin in the 
HSD, which is not observed. 


4.3. Comparison with compact steep-spectrum sources 

In our radio analysis (see Paper II) we And evidence that 
PKS 2004-447 belongs to the class of CSS sources. CSS 
and gigahertz-peaked spectrum (GPS) galaxies are bright ra¬ 
dio sources, which are considered to be young radio galaxies 
based on their peaked radio spectra and compact jet structure 
(e.g., |0’Dea||1998[ |Stangheflini||2003j jPanti et al.jjWTt jRam] 
jdall et al.|2011| and references therein). The same scenario has 
been discussed for y-NLSl due to their low black hole masses 
and high accretion rates (e.g., jPoschini et al.||20T5] l . Besides 
PKS 2004-447, only one other CSS source has been associated 
with a y-ray counterpart in the 3LAC ( jAckermann et al.|2015| l. 
The source thus plays an important role in the study of y-NLS 1 
and CSS sources and a possible connection of these classes. 

The X-ray properties of GPS/CSS quasars have been studied 
by Siemiginowska et al. (|2008| SOS hereafter) and [Kunert- 


Bajraszewska et al. (20141. In most cases, the linear size of 


GPS/CSS sources in radio is smaller than the spatial resolu¬ 
tion of current X-ray instruments. This means that X-ray emis¬ 
sion from the entire complex radio strvicture (i.e., core, jet, and 
hot spots) is contained within the extraction region of the X- 
ray point source. Their X-ray spectra can typically be mod¬ 
eled well by a featureless power law without intrinsic absorp¬ 
tion. The power-law indices range from 1.5 and 2.2, with the 
majority between 1.7 and 2. Their photon index is typically 
close to those of radio-quiet quasars, but there are cases where 
the power law is flat, similar to radio-loud sources. SOS re¬ 
port that for these cases, the X-ray emission can be associated 
with emission from the radio structure. The X-ray absorption- 
corrected luminosities of GPS/CSS sources are in the range 
of (10"^"^ - 10"^®)ergs“'. The X-ray photon index and unab¬ 
sorbed luminosity of PKS 2004-447 thus agree with properties 
of low-powered CSS sources, where the radio jet contributes 
to the X-ray spectrum. Eor some GPS/CSS sources, SOS And 
more complex spectra with evidence of intrinsic absorption of 
~ 10“^'cm“^ and/or a soft excess. The (0.5-2)keV flux con¬ 
tribution of the latter is on the order of 10 '"^ergcm^^ s~'. The 
origin of this emission is not yet well understood. Different the¬ 
oretical predictions include thermal emission through jet-gas in¬ 
teractions with the interstellar medium. The flux of the tentative 
soft excess of PKS 2004-447 in the soft band is consistent with 
findings by SOS. In the context of CSS sources, this may indi¬ 
cate interaction between the young radio jet and the surrounding 
medium, but further studies are necessary. 


4.4. Comparison with other y-NLS1s and biazars 

jPoschini et al.j ( |20I5| ) have studied the multiwavelength prop¬ 
erties of 43 radio-loud NLSl galaxies, including the seven 
y-NLS 1 known so far. PKS 2004-447 is the second clos¬ 
est y-NLSl. It is also the radio-loudest and only southern- 
hemisphere source in the small sample. Besides PKS 2004-447, 
the sources that have been studied in the most detail so far are 
PMNJ094S-H0022, PKS 1502-H036, 1H0323 -h 342, and recently. 
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Fig. 7. X-ray luminosity Lx in the (0.3-10) keV band as a function of 
photon index for y-NLS 1 galaxies in this paper (black filled symbols) 
and the MOJAVE 1 sample (open symbols). 
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Fig. 8. Comparison of (0.3-10) keV X-ray luminosity Lx as a function 
of redshift z for the MOJAVE 1 sample (open symbols) and y-NLS 1 
sources (black filled symbols) discussed in this paper. 


SBS 08464-513 ( 

Abdo et al.||2009b|a| ID’Ammando et al.||2012| 

\2UU\ [Yao et al. 

2015|l. In the radio band, we observe that the 


15 GHz luminosities span two orders of magnitude, and both 
flux and spectral variability are observed in all sources (Paper II), 
though the persistent steep radio spectrum in PKS 2004-447 is 
unique in this sample (Paper II). In this section, we compare our 
results for PKS 2004-447 with characteristics of these sources. 


ture (Zhou et al. 

2007 [Yuan et al. 20081 [Abdo et al.[[2009b 

D’Ammando et a 

[.|2012 2014 2013a 

Paliya et al. 201 3| [Rivers 

et al.|2013l Bhattacharyya et al.|2014 

Yao etal.|2015l. Table [7 


provides results of the simplest best-fit model for each source. 
For comparison, we include results for PKS 2004-447 from this 
analysis averaged over all observations. 

All five y-NLSl objects show moderate flux variability. Only 
the closest, y-NLSl ( IH 03234-342), shows flaring activity on 
timescales of days ( |Paliya et al.||2014{ |Yao et aL]|2015| l. No 
evidence was found of intrinsic absorption in any of these ob¬ 
jects. It is therefore reasonable to conclude that variations in 
flux imply changes in the luminosity and/or spectral shape, 
rather than changes in the obscuration of the primary radiation. 
Spectral variability is observed in two cases (PKS 2004-447, 
IH 0323-1-342) and suspected in two more (PKS 15024-036, 
PMN 109484-0022). 

Results from X-ray analyses reported in the literature mentioned 
above were used to calculate the X-ray luminosities using the 
K correction by, for example, |Ghisellini et al.| ( [2009[ l. Figure [7] 
shows the X-ray luminosities of y-NLSls as a function of the 
photon index. The luminosities span a range of almost two orders 
of magnitude from 10‘*"^ergs“' to lO^^ergs^', while the photon 
indices ranges from 1.2 up to 2.2. However the majority of mea¬ 
surements ofF cluster between 1.5 and 1.8. 

We compare this to the X-ray properties of blazars forming 
a statistically complete sample of radio-selected extragalactic 
jets (MOJAVE 1, [Lister & Homan 20051 that were detected 
by FermifLKT (see Figs. nr and|^. The X-ray data on the 
MOJAVE 1 sources were obtained from single Swift observa- 
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Fig. 9. Comparison of (0.3-10) keV X-ray luminosity Lx as a function 
(0.1-100) GeV luminosity Ly of the MOJAVE 1 sample (open symbols) 
and y-NLSl sources (black filled symbols) discussed in this paper. 


tions ( jChangpOlOj l. X-ray luminosities of the y-detected MO¬ 
JAVE 1 sources span a wider range from lO'^^ergs"' to a few 
10^^ erg s '. In Eigs.|^and[^ we show the X-ray luminosity of 
both samples as a function of the redshift and y-ray luminos¬ 
ity taken from the 2EGL catalog ( [Nolan et al. 112012) . Eor each 
y-NLSl source, we used the weighted mean of the X-ray lu¬ 
minosities from different observations. The blazar distribution 
shows the known sequence from high-luminosity ESRQs de¬ 
tected at high redshifts to lower luminosity BL Lac objects ob¬ 
served at lower redshifts (see, e.g., [Fossati et al.|199^. The X- 
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ray and y-ray luminosities of y-NLSl lie in the same range as 
blazars. More intriguingly, their luminosities fall into a region 
that is occupied by BL Lac objects and low-luminosity FSRQs. 
This result is consistent with the findings of |Paliya et al.| ( |20T3| l, 
who analyzed non-simultaneous SEDs of two y-NLSls, namely 
PKS 2004-447 and PKS 1502-H036. In comparison to the SED 
of a typical BL Lac and ESRQ source, they found that their 
broad-band spectra can be regarded as intermediate between 
both blazar types, although they bear a closer resemblance to 
the SED of ESRQs. Similar results have been obtained by inde¬ 
pendent studies from Foschini et al. (|2015 1 and Sun et al. (2015 i. 


5. Conclusion 

As part of the TANAMI multiwavelength program, we analyzed 
new and archival X-ray observations made by Swift and XMM- 
Newton of the only known southern-hemisphere, radio-loud, y- 
ray emitting narrow line Seyfert 1 galaxy PKS 2004-447 and 
compared the results obtained so far for the four other well- 
studied sources of this type. 

The (0.5-10) keV X-ray spectrum of PKS 2004-447 is best de¬ 
scribed by a power law with Galactic absorption. Its photon in¬ 
dex is consistent with values typically observed for FSRQs. No 
evidence for emission or absorption lines has been found. A pos¬ 
sible soft excess reported in the literature (G06) could not be 
confirmed by our new XMM-Newton observations in 2012, but 
its existence in the data set of 2004 cannot be ruled out. Its flux 
contribution in the soft spectrum of 2004 is less than ~ 20%. 

We observed moderate flux changes with a factor of ~ 3 be¬ 
tween the minimum and maximum flux on timescales down to 
two months. Weak spectral variability was observed on similar 
timescales. However, the hardness intensity diagram shows no 
correlated behavior of count rate and hardness ratio. The linear 
behavior of variations between the (0.5-2) keV and (2-10) keV 
bands indicates that no significant brightness-dependent spectral 
variations occurred. Supported by evidence of a one-sided rela¬ 
tivistic jet in the radio data of PKS 2004-447 (Paper II; see also 
|Orienti et al.||2012) l, we concluded that the observed spectrum 
and variability characteristics can be explained by non-thermal 
emission from the relativistic jet. 

Comparing the X-ray properties of the detected y-ray NLSl 
galaxies, we found that the major difference between these 
sources is their luminosity, which spans a range of almost two 
magnitudes. With the exception of IH 0323-1-342, all spectra 
present a flat and a smooth power-law spectrum. A soft excess 
does not seem to be a frequent spectral component in y-NLSls. 
Their photon indices are clustered between 1.5 and 1.8. 
Continuum variability is a common property of these sources. 
It is most prominent in SBS 0846-1-513, PKS 2004-447, and 
IH 0323-1-342. The analysis of spectral variabilit, however, re¬ 
quires more high-S/N observations. So far, it has only been sig¬ 
nificantly observed in one out of five sources. In comparison 
with typical y-ray detected blazars from the MOJAVE 1 sample, 
we And that X-ray and y-ray luminosities of y-NLS 1 s lie in the 
range of BL Lac objects. 

In a second paper (Paper II), we presented the radio properties of 
PKS 2004-447 based on VLBI and other radio observations. We 
compared the radio spectrum and variability with the same y-ray 
NLS1 galaxies as were presented in this paper. While their X-ray 
spectra show similar steepness and different variability charac¬ 
teristics, the major difference in the radio properties is the radio 
spectrum. Among the five sources, three show a flat radio spec¬ 
trum and one shows a GHz peaked synchrotron (GPS) spectrum. 


while PKS 2004-447 is the only steep spectrum source. The lu¬ 
minosities span a broad range of magnitudes similar to those in 
X-rays. Based on the small linear size of the radio structure and 
the steep radio spectrum, PKS 2004-447 can be considered a 
CSS source, i.e. a young radio galaxy at a moderately small ori¬ 
entation angle to the line of sight. Its X-ray results are in good 
agreement with X-ray properties of this source type and thus 
supports the hypothesis that y-NLSl are young radio sources. 
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Table 7. Overview of X-ray spectral properties of y-ray emitting and radio-loud NLSl galaxies. 


Source 

Redshift 

Model 

F 

^ soft 

F ^ 

A hard 

Eb^ 


Variability® 

Ref 

1H0323-H342 

PKS 1502H-036 

PKS 2004-447 

PMNJ0948H-0022 

SBS0846-H513 

0.061 

0.409 

0.240 

0.585 

0.584 

simple power law 
simple power law 
simple power law 
broken power law 
simple power law 

o 14 + 0.03 
- 0.02 

2.02 + 0.06 

1.7+ 0.2 

1.58 + 0.11 

1 48+004 
°-0.03 

1.47 + 0.25 

1 79+0.09 
^•'^-o.ii 

14.0 + 0.1 

0.3 +0.7 

0.7 + 0.4 

A ^q+0.03 

0.8 + 0.4 

F/S 

F/(S) 

F/S 

F/(S) 

F 

(1) 

(2) 

this paper 

(3) 

(4) 


Notes. Overview of spectral parameters of known y-ray and radio-loud NLSl galaxies for the simplest best-fit models in the literature. Photon 
Index, in the case of PMN J0948-t0022 yi represents the photon index before the break energy. *** Photon Index after the break energy of the 
broken power law. Break energy in units of keV of the broken power law. X-ray flux in units of 10“*^ erg s^’cm^^ in the range of (0.3-10) keV 
in case of PMNI0948-H0022, PKS 1502-H036, and SBS 08464-513, and (0.2-10) keV for 1H0323-H342. <'■> Variability of the X-ray spectrum. Flux 
variability is denoted as F, spectral variability as S. Indications of variability that has not been confirmed are shown in parenthesis. References: 
(l) |Abdo et al.| ( (2009b^ ; (2) |Paliya et al.| ( [2013| l; (3) |D’Ammando et al.| ( |2014^ ; (4) |D’Ammando et al.H2012) 


Kadler M., Ojha R., TANAMI Collaboration 2015, AN 336, 499 
Kalberla P.M.W., Burton W.B., Hartmann D., et al., 2005, A&A 440, 775 
Kara E., Fabian A.C., Cackett E.M., et al., 2013, MNRAS 434, 1129 
Kaspi S., Smith RS., Netzer H., et al., 2000, ApJ 533, 631 
Kellermann K.I., Sramek R., Schmidt M., et al., 1989, AJ 98, 1195 
Komossa S., Voges W., Xu D., et al., 2006, AJ 132, 531 

Kunert-Bajraszewska M., Labiano A., Siemiginowska A., Guainazzi M., 2014, 
MNRAS 437, 3063 

Lister M.L., Homan D.C., 2005, AJ 130, 1389 
Maraschi L., Ghisellini G., Celotti A., 1992, Apjl 397, L5 
Markowitz A.G., Reeves J.N., 2009, ApJ 705, 496 
Marscher, A. P. in Belloni T., 2009, Lect. Notes Phys. 794 
Mason K.O., Breeveld A., Much R., et ah, 2001, A&A365, L36 
Mehdipour M., Branduardi-Raymont G., Kaastra J.S., et ah, 2011, A&A534, 
A39 

Nolan P.L., Abdo A.A., Ackermann M., et ah, 2012, ApJs 199, 31 
O’DeaC.R, 1998, PASP 110, 493 

Ojha R., Kadler M., Bock M., et ah, 2010, A&A p 519, A45 
Orient! M., D’Ammando R, Giroletti M., for the Fermi-LAT Collaboration 2012, 
ai-Xiv astro-ph.HE/1205.0402 

Oshlack A.Y.K.N., Webster R.L., Whiting M.T., 2001, ApJ 558, 578 
Osterbrock D.E., Pogge R.W., 1985, ApJ 297, 166 
Paliya V.S., Sahayanathan S., Parker M.L., et ah, 2014, ApJ 789, 143 
Paliya V.S., Stalin C.S., Shukla A., Sahayanathan S., 2013, ApJ 768, 52 
Piconcelli E., Jimenez-Bailon E., Guainazzi M., et ah, 2004, MNRAS 351, 161 
Randall K.E., Hopkins A.M., Norris R.P, Edwards P.G., 2011, MNRAS 416, 
1135 

Rivers E., Markowitz A., Duro R., Rothschild R., 2012, ApJ 759, 63 

Rivers E., Markowitz A., Rothschild R., 2013, ApJ 772, 114 

Sambruna R.M., Gambill J.K., Maraschi L., et ah, 2004, ApJ 608, 698 

Siemiginowska A., LaMassa S., Aldcroft T.L., et ah, 2008, ApJ 684, 811 

Stanghellini C., 2003, PASA 20, 118 

Striider L., Briel U., Dennerl K., et ah, 2001, A&A 365, L18 

Sun X.N., Zhang J., Lin D.B., et ah, 2015, ApJ798, 43 

Turner M.J.L., Abbey A., Arnaud M., et ah, 2001a, A&A 365, L27 

Turner T.J., Romano P, George I.M., et ah, 2001b, ApJ 561, 131 

Vaughan S., Edelson R., Warwick R.S., Uttley R, 2003, MNRAS 345, 1271 

Vaughan S., Fabian A.C., 2003, MNRAS 341, 496 

Vaughan S., Reeves J., Warwick R., Edelson R., 1999, MNRAS 309, 113 

Vemer D.A., Ferland G.J., Korista K.T., Yakovlev D.G., 1996, ApJ 465, 487 

Wandel A., 2000, New Astr. Rev. 44, 427 

Wilms J., Allen A., McCray R., 2000, ApJ 542, 914 

Yao S., Yuan W, Komossa S., et ah, 2015, ApJ 150, 23 

Yuan W., Zhou H.Y, Komossa S., et ah, 2008, ApJ 685, 801 

Zhou H., Wang T, Yuan W., et ah, 2007, ApJ 658, L13 

Zhou X.L., Zhang S.N., 2010, ApJ 713, LI 1 


Article number, page 12 of 


12 















